/* R. J. Cano, Jan 15 2014 */
polypri(n)=sum(k=1,n,prime(k)*x^(n-k));
vectpri(n)=vector(n,k,prime(k));
vectpriextract(w)=sum(k=1,#w,w[k]*x^(#w-k));
Omega(n)=matrix(n,n,y,x,y>x);
a(n)=(polypri(n)-(vectpriextract(Omega(n)*vectpri(n)~))*(x-1))-n;
